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ABSTRACT 


Four  boundary  value  problems  in  unsteady  simple  shear  flow  are 
considered.  Solutions  for  the  second-order  fluid  approximation  are  listed, 
including  a  solution  which  seems  to  be  new.  For  the  third-order  fluid  approx¬ 
imation,  solutions  are  obtained  for  the  linearized  versions  of  the  various  pro¬ 
blems.  The  Second  Law  of  Thermodynamics,  reduced  to  the  requirement 
that  the  rate  of  deformation  work  be  non-negative,  is  applied  to  these  solu¬ 
tions.  The  signs  of  the  pertinent  second-order  material  constant  is  not 
determined  whereas  the  signs  of  the  pertinent  third-order  constants  are 
determined.  On  the  other  hand,  the  Second  Law  places  restrictions  on  the 
allowable  solutions;  for  some  problems  no  solution  is  valid.  The  implica¬ 
tions  of  these  results  are  discussed. 
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1.  Introduction 


1.1  Simple  fluids  and  n-th  order  approximations 

Noll^  has  defined  a  general  type  of  nonlinear  fluid  with  memory  which  is 
called  a  simple  fluid.  In  this  paper  we  shall  be  concerned  only  with  imcom- 
pressible  fluids  so  that  the  constitutive  equation  takes  the  form 

T  =  -pi  i  (6  w)  (!•« 

s^o 

where  T  is  the  stress  tensor,  p  is  the  indeterminate  pressure,  I  the 

oo  ^ 

identity  tensor,  f  ^  functional  over  the  deformation  history  tensor 

s*?o  ^ 

G  (s).  Throughout  this  paper  we  will  call  T  +  p^the  extra  stress  and  denote 
it  bv  S  .  The  deformation  history  tensor  is  given  by 

whe  re 


and 

F,.  (i-'J  . 

where  ^  is  location  at  time  (t  -  s)  of  the  material  particle  which  is  at  < 
at  time  t.  In  (1.1)  the  functional  must  be  isotropic  in  G(s) 
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W.  Noll,  Arch.  Rati.  Mech.  Anal.  2,  197  (1958). 


and  we  impose  the  conditions  that  if  the  fluid  has  been  at  rest  for  all  time 
the  extra  stress  is  zero,  that  is 


i  (Q)  -  O 


(1.5) 


and  thus  p  is  a  hydrostatic  pressure.  Simple  fluids  include  the  pheno¬ 
menon  of  stress  relaxation  as  well  as  nonlinear  shear  stress  versus  shear 

rate,  and,  unequal  normal  stresses. 

2 

Coleman  and  Noll  have  shown  that  if  the  functional  f  is  n-times  differ- 
entiable  in  a  suitable  way  with  respect  toflie  deformation  history  G  and  if 
G(s)  has  n  derivatives  at  s=  0,  the  following  constitutive  equation  is  a  suitable 
n-th  order  approximation  for  slow  enough  motions: 


T  =  -pi  +  z  i-  r  A-  A-  1 

where  the  summation  is  over  all  sets  of  integers  j..  .  .  j  for  which 

1  X 


(1.6) 


<  i,<  ^  n  +  Jk  ^ 


(1.7) 


The  A  are  Rivlin-Ericksen  tensors  which  have  the  curvilinear -coordinate 


component  form 


rA;), 


iij  Cl)  V  /t  \  '•'p  ■' 

“  (j)  ^ 

O') 


(1.8) 


where,  for  example,  is  the  convariant  derivative  with  respect  to 

the  position  coordinate  x”^  of  the  k-th  component  of  the  i  -th  accelera 


B.  D.  Colman  and  W.  Noll,  Arch.  Rati.  Mech.  Anal.  6,  355  (I960). 
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tion.  The  *,  are  material  coefficients  and  the  relation  (1.6)  is  linear  and 
isotropic  in  each  A^. 

The  precise  manner  in  which  (1.6)  is  an  approximation  to  (1.1)  can  be 
found  in  Ref.  2.  Suffice  it  to  say  here  that  if  the  time  scale  is  slowed  down 
by  a  factor  c<  ^  0  £  ^  ^  I  ,  the  error  made  in  approximating  (1.1)  by 

(1.6)  is  higher  than  order  n  in  oC  .  In  particular  (1.6)  does  not  include 
stress  relaxation  effects;  that  is,  if  the  first  n  time  derivatives  of  the 
motion  vanish  the  stress  reduces  to  the  hydrostatic  pressure. 

A  logical  question  to  ask  is  whether  the  constitutive  equation  (1.6)  can 
be  any  more  than  an  n-th  order  approximation  to  a  simple  fluid.  That  is,  is 
it  reasonable  to  assert  that  (1.6)  could  represent  a  fluid  for  all  histories  no 
matter  what  their  magnitude;  can  we  talk  about  n-th  order  fluids?  This  ques¬ 
tion  will  be  partially  answered  in  the  remainder  of  this  paper. 


1,2  The  Second  Law  of  Thermodynamics 

The  differential  balance  of  energy  can  be  written  in  the  form 


=  tr  (T  d)  - 


(1.9) 


where  ^  rate  of  deformation  tensor,  £  is  the  internal  energy 

per  unit  mass,  q  is  the  heat  flux  vector  and  r  is  the  heat  supply  per  unit 
mass  and  ^  is  the  density.  The  differential  entropy  balance  is  given  by 


f?  -  -2i  -  fi  -  f 


(1.10) 


where  ^  is  the  entropy  per  unit  mass,  &  is  the  absolute  temperature 
and  V  is  the  entropy  production  per  unit  mass.  Following  Coleman  and  Noll' 


B.  D.  Coleman  and  W.  Noll,  Arch.  Rati.  Mech.  Anal.,  13,  167  (1963) 
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we  state  the  Second  Law  of  Thermodynamics  as  the  requirement  that  Y 
be  non-negative  for  all  admissible  thermo-mechanical  processes.  By  ad¬ 
missible  thermo -mechanical  processes  are  meant  processes  which  satisfy 
the  balance  of  linear  momentum  as  well  as  the  balance  of  energy  and  also 
satisfy  all  assumed  constitutive  equations.  The  Second  Law  thus  is  seen  to 
place  restrictions  on  the  allowable  constitutive  equations  and/  or,  as  we 
shall  see,  on  the  allowable  motions.  From  (1.9)  and  (1.10)  we  get  for  the 
Second  Law 


j0y  =  f(e^-i)  ■+  tr(ip)-^2  ye  (lid 

For  an  incompressible  material  (1.11)  reduces  to 

-e)  +  t»"(5 D) d.iz) 

We  can  restrict  our  attention  to  purely  mechanical  variables  by  assum¬ 
ing  that  the  entropy  and  internal  energy  are  constant,  and,  that  the  entropy 
production  due  to  heat  conduction  vanishes,  so  that  (1.12)  reduces  to 

=  tr  (  5D)  >  0  (1. 13) 

which  says  that  the  entropy  production  is  due  entirely  to  the  rate  of  deform¬ 
ation  work.  We  cannot  say  how  much  the  above  assumptions  are  restrictions 
until  we  say  something  about  the  constitutive  equations  for  the  thermal 
variables  ,  A  and  cj  .  Materials  and  conditions  for  which  the  inequality 
(1.13)  holds  are  discussed  by  Coleman  .  We  will  comment  further  on  this  point 


B.  D.  Coleman,  Arch.  Rati.  Mech.  Anal. ,  273  (1962).  The  starting 

point  of  the  relevant  discussion  in  that  refer^ence  is  justified  by  B.  D. 
Coleman,  Arch,  Rati.  Mech.  Anal.,  17,  1  (1964)  which  contains  a  general 
thermodynamic  theory  of  materials  with  memory. 
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in  Section  4.  Such  considerations  are  beyond  the  scope  of  this  paper.  However, 
we  can  note  the  following.  If  there  is  no  heat  flow,  which  is  the  case  for  an 
adiabatic  material,  then  by  (1.9)  and  (1.13)  we  have 

pr  =  -  tr(SD^  <  0 

that  is,  in  order  for  the  balance  of  energy  to  be  maintained  heat  must  be 
withdrawn  at  each  point  at  a  rate  exactly  equal  to  that  determined  b  y  the  rate 
of  deformation  work.  On  the  other  hand,  if  the  entropy  production  due  to 
heat  conduction  vanishes  because  the  fluid  is  such  a  good  conductor  that  V 0 
is  zero,  then  by  (1.9),  if  f  =  o,  we  require  that 

V-q  -  tr(rj) 

I  '  /V  ^ 

that  is,  q  is  determined  by  the  rate  of  deformation  work.  (Note  that  in  the 
latter  case  q  is  indeterminate  in  the  sense  that  it  is  not  determined  from  a 
constitutive  equation.  The  situation  is  quite  analogous  to  the  indeterminate 
pressure  in  incompressible  fluids). 

1,3  Unsteady  simple  shear  flows:  boundary  value  problems 

An  unsteady  simple  shear  flow  is  a  flow  whose  velocity  field  V  can  be 
represented  in  a  Cartesian  rectilinear  coordinate  system  by  the  following 
components 

U-U(y^'t)  V'O  Yt  -  O  (1.16) 

The  differential  equation  form  of  the  balance  of  linear  momentum  can  be 
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expressed  in  vector  notation  by 


IZ  ■>  fk  ■  ('•"> 

where  V  is  the  derivative  operator  and  L  is  the  external  body  force 
per  unit  mass.  We  will  assume  that  can  be  derived  from  a  potential 

X  by  means  of  f^  =  -VX,'Vfe  then  introduce  in  the  usual  way  the  modi- 
fied  pressure  (f)  by  means  of 


^  ^  -f-  X  (1.18) 

Then,  in  terms  of  the  extra  stress  5  the  balance  of  linear  momentum 
becomes 

§,Y  - -  fi  (1-19) 

Since  the  velocity  is  assumed  to  be  independent  of  X  and  x  it 
follows  that  the  Rivlin-Ericksen  tensors  are  also  independent  of  x  and  z 
and  also  along  with  (1. 16)  it  follows  that 

5  xz  "■  -  5zz  ~  ^  (1.20) 

It  then  follows  that  the  components  of  (1. 19)  reduce  to 

~  ^  ~  tA 

dy  5y^  -  ^yCi>  =  0 

d^4)  *  o 
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(1.21) 

(1.22) 

(1.23) 


We  see  that  the  modified  pressure  must  be  independent  of  z.  and  from  (1.21) 
and  (1.22)  that 


0  =  5yy  ^  ^  U)  (1*24) 

and  hence  (1.21)  reduces  to 

ay  5;,^  =  w  ip  it)  (1.25) 

where  ip  (t)  is  the  modified  pressure  gradient  in  the  x-direction.  The 

analysis  of  this  paragraph  is  a  generalization  of  one  contained  in  Markovitz 
5 

and  Coleman. 

We  now  list  several  boundary  value  problems  in  unsteady  simple  shear 

flow: 

Problem  1;  We  consider  the  boundary  conditions 

u(o,t)  •  U  u(L,l)-0  (1.26) 

along  with  the  assumption  that  the  modified  pressure  (j)  has  no  gradient 
in  the  x  -  direction  so  that  by  (1.24) 

0(t)=O 

For  arbitrary  initial  conditions  this  problem  is  transient  as  well  as  unsteady. 
^  H.  Markovitz  and  B .  D.  Coleman,  Phys.  Fluids,^,  833  (1964). 
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For  this  problem  we  make  the  following  non-dimensionalizati' 


U  = 


t  =  -id 
fL- 


y't 


(1.28) 


and  the  boundary  conditions  (1,26)  reduce  to 


=  i  n(],0-o 


(1-29) 


The  rest  of  the  problems  we  consider  are  unsteady  but  oscillatory  and 
so  do  not  require  initial  conditions. 


Problem  2: 


=  UcosSlt  u(Lj)  = 


(1.30) 


This  time  we  make  the  following  non-dimensionalization 


-  u 
U  - 


t  =  -nt 


^  ^  L 


which  reduces  (1.30)  to 


(1.31) 


(  0 , t  )  =  COS  t 


(l,  ^  )  =  0 


(1.32) 


Problem  3:  Consider  a  variation  of  the  preceding  problem  in  which 
L  -  c?o  >  that  is  we  have  the  boundary  conditions 


(o.  t)  =  Uc05  i?  t  u(o<},t)  =  0 


(1.33) 
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We  introduce  the  non-dimensionalization 


t  =  jit  y^/J^ 

which  gives  the  dimensionless  boundary  conditions 

5  (0,n  =  cost  u(^,i)-  0  (1.35) 

Problem  4;  A  final  problem  in  unsteady  simple  shear  flow  has  the 
boundary  conditions 

ua.tl  -  0  (1.36) 

with  the  pressure  gradient  in  the  x-direction  given  by 

=  <4  COi-Ot  (I*  37) 

With  the  non-dimensionalization 

5-f"^  i-^t  y-t  (?-!■.  (1.38) 

(1.36)  and  (1.37)  reduces  to 

u(l,i)=  a  (-1,  i)  =  0  (1.39) 

and 

p a)  *  I 
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2.1  Solutions 


For  n*;?,  equation  (1.6)  reduces  to 

J  =  "pI  (2-1) 

where  the  fact  has  been  used  that  for  incompressible  fluids  tr  Ai  "  0 
and  also  all  isotropic  tensors  have  been  absorbed  in  the  indeterminate 
pressure  term  .  For  n  =  l  we  would  have  only  the  linear  term  yuA| 

which  comprises  the  classical  linear  theory  of  viscous  fluids.  Each  of  the 
second-order  terms  is  nonlinear:  the  first  because  A.  is  squared  and  the 
second  because  contains  a  product  term  as  given  by  (1.8).  The  co¬ 
efficient  JA  is  the  usual  linear  viscosity  and  and  are  second- 

order  material  coefficients. 

For  unsteady  simple  shear  flow  Coleman  and  Noll  in  Ref.  2  derived  the 
following  non- zero  components  of  the  extra  stress: 

3yy  “  (/^/  ^  f^t  ^  ^y 

and  from  (1.25)  the  governing  differential  equation  for  the  velocity  profile: 

y.  Uyy  +  jU^^l^yyt  =  ^  (2.3) 

where  the  subscripts  on  U  indicate  partial  differentiation. 

It  was  noted  by  Coleman  and  Noll  that,  rather  surprisingly  for  a  non¬ 
linear  theory,  the  velocity  profile  is  governed  by  a  linear  equation.  However 
from  (2.2)  it  is  seen  that  the  normal  stresses  are  nonlinear  as  well  as  unequal. 


We  now  present  the  solutions,  for  the  second-order  approximation  theory, 
to  the  boundary  value  problems  listed  in  Section  1.3. 

Problem  1:  Using  (1.27)  and  (1.28)  equation  (2.3)  reduces  to 


^yy 


.  <*)_ 

A  upyt 


where 


is  the  dimensionless  group 


The  steady- state  solution  of  (2.4)  is 


(2.4) 


(2.5) 


a  I  -  y  (2.6) 
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and  the  unsteady  solution  which  can  be  found  by  separation  of  variables  , 
after  nondimensionalizing  in  the  above  way,  is 


n-  I 


(2.7) 


We  see  that  a  departure  in  the  initial  velocity  distribution  from  the  steady 
state  is  represented  by  a  Fourier  series.  On  examination  of  (2.7)  we  note 
that  for  positive  that  is  for  a  positive  value  of  the  material  constant 

I  every  term  of  the  series  is  convergent,  whereas  for  negative 
only  terms  for  which 


y\  < 


(2.8) 


T.  W.  Ting,  Arch.  RatL,  Mech.  Anal.  14,  1  (1963) 
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are  convergent. 

Using  (1,27)  and  (1.31)  equation  (2,3)  reduces  to 


^  //  ^  ^  yy  ^  ~ 


(2.9) 


u) 


We  see  in  this  case  that  there  are  two  dimensionless  groups:  /\  as  given 
by  (2 . 5)  and  /\  by 


A  (2.10) 

8 

The  solution,  obtained  by  Markovitz  and  Coleman  ,  after  non-dimension- 
alization  in  the  above  way  is 

U  =  ^e.  [  U{y)c'^  j  (2.11) 

where 

u(y)  =  (2.12) 

s/hJ? 


and 


m 


[  i/(A  -  (  A” )] 


(2.13) 


Problem  3:  Using  (1.27)  and  (1.34)  equation  (2.3)  takes  the  form 


7 


8 


See  also  B.  D.  Coleman,  R.  J.  Duffin  and  V.  J.  Mizel,  'Instability 

^xtx 


and  Uniqueness  Theorems  for  the  Equation  u  =  u  -  u  on  a  Strin" 

_ •_  .A„  ,  r,  ,  .  ,  t  XX  ’ 


to  appear  in -Arch.  Rati.  Mech.  Anal. 


H.  Markovitz  and  B.  D.  Coleman,  Adv.  Appl.  Mech.,  Vol.  8,  Academic 
Press,  N.  Y  .  (1964). 
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(2.14) 


Uyy  4-  A  ‘^yyf  ~ 


a‘  * 

where  now  the  dimensionless  group  /\  is  given  by 


(2.15) 


The  solution  to  this  problem  was  obtained  by  Markovitz  and  Coleman  in 
Ref.  5  and,  after  the  above  non-dimensionalization,  is 


uCyj)  '  <f  cos  {I  - 


(2.16) 


where 


A.  B  = 


Z[)  -t 


(2.17) 


The  parameters  A  and  B  are  always  real  and  positive  regardless  of  the 
sign  or  magnitude  of  /\ 

Problem  4:  Using  (1.38)  equation  (2.3)  reduces  to 


A  ^yy  +  A^^^yyi  “  Uj  -f  I  -/■  (fjcost 


(2 . 18) 


where  A  is  given  by  (2.10)  and  A  by  (2.5).  The  non-transient  solu¬ 
tion  to  this  problem  can  be  found  by  assuming  a  solution  of  the  form 


U(y.l)  = 


(y)  ■+  (?)  <?'  ^ \ 


(2.19) 
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One  then  finds  that  the  solution  is  given  by 


‘A;  (7“-') 


«/ (y^  =  ‘&)l'  ~ 


cosU  jfty 
cosh  tn 


(2.20) 


(2.21) 


where  is  given  by  (2.13). 


2.2  Application  of  the  Second  Law 


Substituting  the  constitutive  equation  (2.1)  of  the  second-order  approxima 
tion  into  the  Second  Law  inequality  (1.13)  we  find  that 

^eY  ^  juir  A,  /^TtrA,  ti- {A,  A,)  >  q  ^2.22) 

A  solution  to  Sec.  2J  which  satisfies  the  thermal  assumptions  leading  (1.13) 
must  also  satisfy  the  inequality  (2.22). 

For  unsteady  simple  shear  (2.22)  reduces  to 


^  2/AUy  -+  > 


(2.23) 


We  see  from  (2.23)  that  for  steady  motion 


/<  >  O  (2.24) 

which  is  the  well-known  result  that  the  viscosity  of  the  classical  linear 
theory  must  be  positive. 

Consider  a  flow  in  which  the  velocity  gradient  oty  =  O  at  some  instant 
of  time  r  for  some  y  .  Then  is  a  minimum  at  that  instant  and 
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0  and  by  (2.23)  =  0  at  such  an  instant.  But  also  by  (2.23)  5)6/ 

must  be  a  minimum  when  ^Q^-0  so  that  ■  0  ,  that  is 


A  If-  '  o 


(2.25) 


There  are  two  possibilities:  =  O  ,  ^  (iAy)/di  ~  0 

which  is  equivalent  to  ^Uy/^t  ~  O  •  Now  assuming  that  is 

unequal  to  zero  in  order  that  we  are  dealing  with  the  second-order  approxi- 
mational  rather  than  the  linear  approximation,  we  have  the  conclusion  that 
a  flow  for  which  and  ^Uy/bt  ^  O  violates  the  Second  Law  of 

Thermodynamics  in  the  form  (1.13).  In  the  light  of  this  remark  we  now 
examine  the  four  solutions  for  unsteady  simple  shear  flow  contained  in 
Section  2.1. 

Problem  1  :  The  velocity  gradient  calculated  from  (2.7)  is 


Uy  =  -/ 


Let  us  suppose  that  there  is  ay#  at  a  certain  instant  of  time  aich 


U  5  (  Xo  ^  =  0  =  “*  “•'2 


(2.27) 


where 


>'ir  A„  e.f  wy. 


(2.28) 


that  is,  the  velocity  profile  has  a  maocimum  or  minimum  at  a  certain  instant 
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of  time.  Now  also  at  and  we  have 

/_J2lZL_^B.  (2.29) 

V  H- n^TT^ 

l  —  *  o  n  —  I 

It  can  be  shown  that  if  (2.27)  holds  then  (2.29)  does  not  vanish.  Thus,  for 
this  problem  flows  for  which  the  velocity  gradient  vanishes  are  not  allowed. 

It  IS  clear  that  if  /\  is  positive  and  the  initial  profile  is  monotomic  then 
zero  velocity  gradients  will  never  develop.  On  the  other,  it  is  clear  that 
if  A  is  negative  and  terms  for  which  (2.8)  is  violated  are  admitted,  then 
such  terms  will  ultimately  lead  to  zero  velocity  gradients.  Therefore  we 
conclude  that  for  negative  A  \  terms  for  which  (2.8)  is  violated  not  only 
give  rise  to  divergent  flows  but  also  violate  the  Second  Law  of  Thermo¬ 
dynamics. 

More  information  about  restrictions  on  the  motion  can  be  found  by 
applying  the  Second  Law  to  a  flow  for  which  the  summation  in  (2.7)  contains 
only  one  term,  that  is 

S  >  l-y  +  A,  gxp  [-  ('y- Wj;  (2.30) 

Substituting  in  the  appropriate  dimensionless  form  of  (2.23)  we  get 

-+  i  >  0  (2.31) 


where 


a  .  0  ,1  A„  (2.32) 


A 


(0 


1  -+  A 


«) 


(2.33) 
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Now  the  left-hand  side  of  (2.31)  is  a  quadratic  form  in  a  and  it  has  the 
discriminant  which  is  positive  or  zero.  For  0  then  the  inequality 
is  violated  for  some  values  of  a  whereas  for  b  =  0,  that  is  0  ,  the 

inequality  is  satisfied.  The  question  then  is:  for  a  given  .  for  what 

values  of  n  and  An  inequality  satisfied  for  all  y  ?  For  /\  ^ 

positive  we  have 

0<  b  <  I  for  A  >  0  (2.34) 

which  meains  that  the  entropy  production  is  negative  when  a  lies  in  the 
range 


- L-  <  a  <  -I 

I  -  b 


(2.35) 


that  is 


-  (I  -t-  n^TT^  <  d  <  -  > 

From  (2.32)  and  (2.36)  we  see  that  for 

hit  (  AA  ^  I 


(2.36) 


(2.37) 


a 


does  not  enter  the  range  of  (2.36).  Another  way  of  stating  (2.37)  is  that 


in  the  form 


Next  let  us  consider  negative  A  .  We  then  write  b 


I  -  n*  «  H 


(2.39) 


Now  for  n  restricted  by  (2.8)  b  is  negative  and  hence  (1-b)  is  positive  so 
that  again  the  entropy  production  is  again  negative  over  a  finite  range  of  a, 


this  time 


-I  <  a  <  -(i  -  n*  TT^  I  ) 


(2.40) 


From  (Z,32)  and  (2.40)  we  see  we  must  have  the  condition 


MTT  I  I  £  I  -  h'-  11  I  I 


(2.41) 


which  is  a  more  severe  restriction  than  (2.38).  This  can  also  be  put  in 


the  form 


-Z  4  I  <  Uy  <  -  A^^^l 


(2.42) 


Problem  2;  F rom  (2 . 11)  we  have 


(2.43) 


where  from  (2.12) 


Uy  - 


ry\  cosh  w  ( i-y) 


^ihh  m 


(2.44) 
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Recalling  that  W  as  given  by  (2.13)  is  complex  we  see  that  there  are  many 
possibilities  for  to  be  zero.  It  is  clear  from  the  boundary  conditions 

of  the  problem  that  for  a  given  y  the  velocity  gradient  is  periodic  with  the 
same  frequency  as  the  oscillating  wall^  Now  from  (2.43)  we  have 

II?  =  gl|i  Wy  I  (2.45) 

Now  let  and  therefore 

up  =  +  (2.46) 

When  UpsO  we  have  D-n/Z  modulo  IT  .  Now 

\  ^  )  (2.47) 

which  is  not  zero  when  ?  )  ■*  modulo  Tl"  ,  Therefore  this  problem 

violates  the  Second  Law  for  all  values  of  the  parameters  of  the  problem. 

Problem  3:  From  the  solution  (2.16)  and  (2.17)  we  get  that  the 

velocity  gradient  vanishes  for  all  ^  which  satisfy 

taH  (T,-  By.).  (2.48) 


For  such  (  Yo  ,  ^0^ 


'^y  ^  cos  (t^- 

at  >5^ 
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Now  from  (2.48)  we  see  that 


“  B  ^  tnodolo  IT 

so  that  (2.49)  would  never  vanish.  Thus  the  solution  for  this  problem  also 
violates  the  Second  Law  for  all  values  of  the  problem  parameters. 

Problem  4;  From  (2.19),  (2.20)  and  (2.21)  we  have 

Wy  =  -b  I  Uj^  ^  (2.50) 

and 


~  (2.51) 

where 


yA 

u,^  = 


-  t  —  ^ 

Cosh  m 


(2.53) 


At  y  =  0  both  Up  and  'buy/ht  vanish  so  that  at  least  on  the  center  line 
the  Second  Law  is  not  violated.  However,  there  may  be  other  y  's  for 
which  the  velocity  gradient  vanishes,  but  bvip/^T  does  not.  The 
answer  to  that  question  appears  to  be  a  complicated  algebraic  problem,  but 
the  author  suspects  that  this  flow  will  be  found  not  to  violate  the  Second  Law 
of  Thermodynamics. 
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3.  Third-order  approximation 


3.1  Solutions  to  linearized  problems 
For  n  =  3  equation  (1.  6)  yields 

^  -  -pi  -t^Ai  •+/^/A,  A^-/- 

-+  )^'  /^i  (3.1) 

+  ("A.  6'  3  A? A, 

where  in  addition  to  the  linear  and  second-order  terms  of  (2.1)  we  have  five 
new  terms,  each  with  a  new  material  constant. 

For  unstesidy  simple  shear  flow  the  extra  stress  components  are 

Sxy  =  fA.  Uy  -f-  JUi^^Uyi  -i-  ^yit 


Sxz  =  Syz  -  O 

Sxx  -  ^y  -i-  2^^^  Uy  ^yt 

syy  =  (/.?'  t  2 /“J  “r 

Szz  =  0 


(3.2) 


whe  re 


fd 


yAj 


Substituting  the  first  equation  of  (3.2)  into  (1.21)  we  get 


(3.3) 


/^z  ^yy't 


ty)  Ulyy 


(3) 

V-  /^r  ^yy  ^ 


ip(t) 


(3.4) 
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We  note  immediately  that  this  governing  differential  equation  is  nonlinear 
and  also  that  we  now  have  second-order  time  differentiation.  We  present 
here  solutions  to  linearized  versions  of  the  boundary  value  problems  listed 
in  Section  13. 

Problem  1:  Using  (1.27)  and  (1^28)  equation  (3.4)  reduces  to 


^yy  ■+  A  +3  A,  (U^fu^y  (3 

where  in  addition  to  A  given  by  (2. 5)  we  have  two  more  dimensionless 
numbers 


(3.6) 


Again  we  have  the  steady -state  solution  (1  -  y  ).  Since  (3.  5)  is  non¬ 
linear  we  are  not  surprised  to  find  that  an  exact  solution  by  analysis  is  not 
readily  obtained.  For  example,  the  method  of  separation  of  variables  which 
can  be  used  to  find  the  solutionto  the  linear  equation  (2.4)  does  not  work  in 
this  case.  We  therefore  settle  for  a  perturbation  analysis  of  the  steady  state 
solution;  we  assume  a  series  solution  of  the  form 


u  =  i-x  +  (3.7) 

where  e  is  a  parameter  which  tends  to  zero.  Upon  substitution  into 
(3.  5)  and  equating  coefficients  of  6-  we  have 


('r+jA‘;’:)v  4  4  A 


2  A/*'  '  ''i 


(3.8) 
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which  is  a  linear  differential  equation  for  the  first  term  V  Ty,  t  )  in  the 
series  of  (3.7).  The  velocity  perturbation  V  must  satisfy  the  boundary 
conditions  that  it  vanish  on  both  boundaries. 

By  separation  of  variables  we  obtain  the  following  solution  of  (3.5): 


Because  of  the  second-time  derivative  we  must  specify  the  acceleration  as 
well  as  the  velocity  at  time  t  =  0  in  order  to  determine  a  solution. 

A  (2)  (^3) 

A  and  /\  is  the  above 

j  '  'I  ^ 

solution  convergent?  Off  hand  it  does  not  appear  that  the  condition  (2.8) 
for  negative  A*'  automatically  carries  over  to  the  stability  analysis  of 
the  third-order  theory.  The  conditions  for  which  the  solution  (3.9)  and 
(3.10)  converges  or  diverges  are  summarized  in  Table  1.  We  note  that  for 
stability  ■+  A”0  ^  and  (l-f  )  must  all  have  the  same 

sign.  If  is  negative  then,  by  (3.6)  ((  +  3  A^*M  could  change  sign 

as  U/  L  is  varied.  However  and  4  /^> )  cannot  change  sign 

with  U  /Li.  Therefore  in 
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order  to  have  stability  for  negative  the  quantity  U/  L  is  restricted 

(3) 

to  ranges  for  which  (1  +  3  A,  )  does  not  change  sign.  Such  ranges  are 


o<U  < 


for  1+3 A.  >0  (3.11) 


3 


57^ 


(3.12) 


Problem  2:  Using  (1.27  and  (1.31)  equation  (3.4)  reduces  to 


A  ■+A\^UT,yT  ■+ SA,  (Uy) 


(3.13) 


In  addition  to  the  dimensionless  groups  /\  and  /\  given  by  (2, 10)  and 


(2.5)  respectively  we  have  the  two  new  groups 

A  '»  i;'"  /u 


'  f>L‘JZ(L 


A  -- 


(3.14) 


For  Hf  =  0  ,  and  hence  /\,  =  O  ,  there  is  no  motion,  after  any  transient 
motion  has  died  out.  We  therefore  consider  small  motions  by  assuming  the 

A  (3) 

following  series  in  /\|  : 


(y,f;  =  A,  V +  o[('A;’'y] 


(3.15) 


On  substituting  (3.15)  into  (3.13)  and  equations  coefficients  of  /\,  we  get 
the  following  linear  equation  for  V 


A  ypy 


(3.16) 
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The  solution  of  (3.16)  is  just  a  slight  modification  of  that  for  the  corres¬ 
ponding  solution  for  the  second-order  approximation  theory ;  v  is  given 
by  the  right-hand  side  of  (2,11)  and  (2.12)  but  (2.13)  is  replaced  by 


(3.17) 


Problem  3.  Using  (1«27)  and  (1.34)  equation  (3.4)  reduces  to 


(2)  .  <2^  ^ 

yy  +A  dypt  -+-  3  A,  (  Up)  (3.18) 


2.'^  y/ 


k(Z)  a  (3)  (?) 

where  /\  is  given  by  (2. 15)  and  /\  and  /\  by 


(-3)  _(?) 


A  = 


(3.19) 


As  in  the  previous  problem  we  look  for  a  solution  f  or  small  motions  by  means 
of  (3.15).  The  resultant  linear  equation  for  is 


Ay  ^  A  Vyyi  A,  VppR  -  v^- 


(3.20) 


The  solution  to  this  problem  is  given  by  the  right-hand  side  of  (2.16)  with 
(2.17)  replaced  by 


A,  8  - 


(3.21) 


The  parameters  A  and  B  are  always  real  and  positive  regardless  of  the 
signs  or  magnitudes  of  and 
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Problem  4:  Using  (1.38)  equation  (3.4)  reduces  to 


+3A'*(5P"%  =  “f  (3-221 

where  A  is  given  by  (2. 10)  and  /\  by  (2.5)  and  /\,  and  /\  are 


(3) 


given  by 


(3)  -(?),/.^ 


la;  —  vA/, 

A 

f 


'  -  VL'' 
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L 


(3.23) 


Again  we  look  for  a  solution  for  small  motions  by  means  of  (3.15).  The 
resultant  linear  equation  for  V  (y is 


A  v-v.-  +  X‘^v„f  -<■  Ki'^yyV  =  Vf  -t  I 


y/ 


(3.24) 


The  solution  to  this  problem  is  given  by  the  right-hand  side  of  <2.19) 
by  (2.20)  and  (2.21)  but  (2.13)  is  replaced  by  (3.17). 

3.2  Application  of  the  Second  Law 

Substituting  the  constitutive  equation  (3.1)  of  the  third-order  approxima¬ 
tion  into  the  Second  Law  inequality  (1.13)  we  find  that 

+  y^rtrA^  H-  H- 

-hya!^\ir  )  -+  /<j.  it~Al  ■+  (3.25) 

^  2/-"'  (a'  aX  +  A-'  ~ 
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For  unsteady  simple  shear  flow  (3.2  5)  reduces  to 

Z  ■+  UyUytf/A  Uy  UyUyU  Z  0  (3.26) 

We  see  from  (3.26)  that  for  steady  motion  we  must  have 


>  0 


(3.27) 


as  well  as  (2.24). 

As  in  Sec.  2.2  we  again  find  by  examining  (3.26)  that  the  rate  of 
change  of  the  velocity  gradient,  that  is  hUy/dt  ,  must  be  zero  when  the  velo¬ 
city  gradient  Uy  is  zero  in  order  that  (3.26)  not  be  violated.  We  now 
investigate  what  this  and  (3,27)  mean  in  terms  of  the  four  solutions  of  the 
previous  section. 

A 

Problem  1;  By  (3.6)  equation  (3.27)  implies  that  A,  must  be  non¬ 

negative,  that  is 

.  (3) 

A,  -  0  (3.28) 

Therefore  by  Table  1,  in  order  to  have  stability  (2*8)  must  be  satisfied 
and 

A2  ^  O  (3.29) 

That  is  we  have  the  same  condition  on  ^  and  the  wave  number  n  as  we 
had  for  the  second-order  approximation  and  also  the  requirement  that 
be  positive  for  stability.  It  is  interesting  to  note  from  (3.10)  that  perturba¬ 
tions  for  which 


4Ar(,^3Ar)>fe-A“) 


(3.30) 
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decay  sinusoidally  rather  than  purely  exponentially. 

It  likely  can  be  shown  that  as  in  the  case  of  the  second-order  approx¬ 
imation  for  this  problem  the  velocity  profile  for  a  perturbed  steady  flow 
must  be  monotonic  in  order  that  the  Second  Law  is  not  violated.  It  likely 
can  also  be  shown  that  violation  of  the  stability  conditions  also  violates  the 
Second  Law. 

Problems  2,  3,  and  4:  Since  the  solutions  of  the  linearized  versions 

of  these  problems  differ  from  the  solutions  of  the  corresponding  second-order 
problems  by  only  unimportant  changes  in  constants  and  since  the  Second  Law 
criterion  regarding  zero  velocity  gradients  is  unchanged,  all  of  our  results  of 
Sec.  2.2  for  these  problems  carry  over  here.  That  is,  the  solutions  for 
Problems  2  and  3  are  not  allowed  for  the  third-order  approximation  whereas, 
as  far  as  we  have  investigated,  the  solution  of  Problem  4  is  allowed. 

4.  Discussion  and  conclusions 

We  have  found  that  the  Second  Law  of  Thermodynamics  in  the  form  of 
the  requirement  that  the  rate  of  deformation  work  be  non-negative  places 
restrictions  on  allowable  motions  for  the  second  and  third-order  approxima¬ 
tion  to  fluids  as  well  as  placing  restrictions  on  the  material  constants  appear¬ 
ing  in  these  theories. 

Let  us  first  discuss  the  restrictions  on  material  constants.  In  addition 
to  the  classical  result  that  the  linear  material  coefficient  jU.  is  positive, 
we  found  from  the  Second  Law  that  the  third  order  material  constants 
and  must  be  positive.  On  the  other  hand,  we  see  that  the  sign  of 

the  second-order  material  constant  which  appeared  in  the  analysis  for 

unsteady  simple  shear  flow,  is  not  determined.  It  may  be  possible  to  find 
restrictions  on  the  signs  of  the  second-order  coefficients  by  examining  flows 
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other  than  unsteady  shear.  It  is  interesting  to  note  that  other  evidence 
indicates  that  is  non-zero  and  negative.  Coleman  and  Markovitz^ 

have  shown  that  if  one  assumes  on  the  'basis  of  thermodynamic  intuition" 
that  the  stress  relaxation  function  of  linear  viscoelasticity  is  positive  for 
all  times,  then  is  negative.  Furthermore,  experimental  determina¬ 
tion  of  ,  in  particular  by  Markovitz  and  Brown^^,  have  yielded  only 
negative  values  thus  far. 

We  turn  now  to  the  restrictions  placed  on  allowable  motions.  It 
turned  out  that  for  the  oscillating  wall  problems,  Problems  2  and  3, 
the  solutions  for  the  second  and  third-order  approximation  theories^ were 
not  valid  for  any  values  of  the  parameters  of  these  problems.  Problem  1, 
the  problem  where  one  wall  is  moving  with  constant  velocity,  only  has  valid 
solutions  for  restricted  initial  conditions:  the  velocity  profile  must  be 
monotonic  and,  if  is  negative,  cannot  contain  harmonics  above  a 

certain  value  which  is  dependent  on  | .  On  the  other  hand,  the  solu¬ 
tion  to  Problem  4  when  the  flow  is  driven  by  a  pulsating  pressure  gradient, 
is  valid  as  far  as  we  have  checked.  It  is  interesting  to  note  in  Problem  1 
that  the  conditions  for  stability  are  consistent  with  satisfying  the  Second 
Law.  From  both  stability  and  thermodynamic  considerations,  we  conclude 
that  the  second  and  third-order  theories  are  just  approximation  theories  and 
cannot  be  used  indiscriminately  as  theories  for  any  physically  well-posed 
problem.  If  such  things  as  second-and  third-order  fluids  existed  in  their 
own  right  we  would  expect  to  get  valid  solutions  for  Problems  2  and  3,  and. 
Problem  J.  for  all  initial  conditions.  On  the  other  hand  the  linear  theory  can 
be  used  as  a  fluid  theory  in  that  solutions  for  any  physically  well-posed  problem 
do  not  violate  the  Second  Law.  In  1951  Truesdell^^  suggested  that  the  inequality 

9 

B.  D.  Coleman  and  H.  Markovitz,  J.  App.  Phys.,  35,  1  (1964). 

H.  Markovitz  and  D.  R.  Brown,  Trans,  Soc.  Rheol. ,  ^,  137  (1963). 

11 

C.  Truesdell,  J.  Math.  Pure  Appl. ,  (9),  30,  111  (1951). 
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(1.13)  be  a  restriction  not  on  constitutive  equations  but  rather  on  allowable 
motions,  Coleman  in  Ref.  4  advanced  the  point  of  view  that  (1.13)  must  hold 
for  certain  motions,  namely  those  for  which  the  reduction  from  (1.12)  to  (1.13) 
is  valii;  if  (1.13)  be  violated  in  one  of  these  motions,  then  the  constitutive 
equation  should  be  rejected..  At  least,  the  constitutive  equation  should  be 
rejected  for  those  motions  for  which  the  Second  Law  is  violated. 

It  may  be  that  the  Second  Law  inequality  (1.13)  as  a  requirement  on 
unsteady  problems  considered  in  this  paper  fqr  the  second-and-third-order 
fluid  approximations  is  inconsistent  with  the  simultaneous  approsd  mations 
appropriate  to  the  thermal  variables.  It 'may  be  necessary  to  set  to  zero 
certain  second  and  third -order  thermal  material  constants  in  order  to 
reduce  to  (1.13).  It  is  intended  to  investigate  this,  point  in  the  near  future. 

A  flow  is  said  to  be  a  helical  flow  if  there  exists  a  cylindrical  coordinate 
system  (r,  G  z)  which  the  physical  components  of  the  velocity  have  the 
form 

In  Ref.  5  it  is  developed  that,  for  the  second-order  approximation,  the  com¬ 
ponents  W  and  u  of  the  velocity  field  are  determined  by  two  separated, 
linear,  third-order  partial  differential  equations.  Problems  corresponding 
to  Problems  2  and  4  are  then  solved  for  in  the  case  of  Poiseuille  flow  (w*o) 
and  corresponding  to  Problem  ^  in  the  case  of  Couette  flow  =  o)  •  Explicit 
solutions  are  obtained  when  appropriate  "small  gap"  approximations  are  made. 
The  Second  Law  takes  the  following  form  for  helical  flow  of  the  second-order 
fluid  approximation: 

pe  Jf*  =  -t  ul]  -h  ul]  ^  O  (4.2) 
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For  Poiseuille  flows  (4.2)  reduces  to 


pdV  =  -+  >  0  (4.3) 

and  for  Couette  flows  (4.2)  reduces  to 

^dy  =  ^  (4.4) 

It  would  appear  that  the  results  we  have  obtained  from  applying  the  Second 
Law  to  unsteady  simple  shear  flows  will  carry  over  to  the  analogous  pro¬ 
blems  in  Poiseuille  and  Couette  flowsof  Ref.  5. 
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